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Abstract 


A search for the decays of the Higgs and Z bosons to J/V 7 and T(nS )7 (n = 1,2,3) is performed 
with pp collision data samples corresponding to integrated luminosities of up to 20.3 fb^ 1 collected at 
Vs = 8 TeV with the ATLAS detector at the CERN Large Hadron Collider. No significant excess of 
events is observed above expected backgrounds and 95% CL upper limits are placed on the branching 
fractions. In the J/Vt final state the limits are 1.5 x 10 -3 and 2.6 x 10 -6 for the Higgs and Z boson 
decays, respectively, while in the T(1S, 25,35) 7 final states the limits are (1.3,1.9,1.3) x 10 ~ 3 and 
(3.4,6.5,5.4) x 10~ 6 , respectively. 
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A search for the decays of the Higgs and Z bosons to J/i/j'y and T (nS) 7 (n = 1, 2,3) is performed 
with pp collision data samples corresponding to integrated luminosities of up to 20.3 fb ~ 3 collected at 
yfs = 8 TeV with the ATLAS detector at the CERN Large Hadron Collider. No significant excess of 
events is observed above expected backgrounds and 95% CL upper limits are placed on the branching 
fractions. In the J/ip 7 final state the limits are 1.5 x 10 -3 and 2.6 x 1CP 6 for the Higgs and Z boson 
decays, respectively, while in the 2 S, 3 S) 7 final states the limits are (1.3,1.9,1.3) x 10 -3 and 

(3.4, 6.5, 5.4) x 10 —6 , respectively. 


Rare decays of the recently discovered Higgs bo¬ 
son (BE] to a quarkonium state and a photon may offer 
unique sensitivity to both the magnitude and sign of the 
Yukawa couplings of the Higgs boson to quarks 0H3- 
These couplings are challenging to access in hadron col¬ 
liders through the direct H —► qq decays, owing to the 
overwhelming QCD background [7]. 

Among the channels proposed as probes of the light 
quark Yukawa couplings 01 m, those with the heavy 
quarkonia J/ij) or T (nS) (n = 1,2,3), collectively de¬ 
noted as Q , in the final state are the most readily ac¬ 
cessible, without requirements for dedicated triggers and 
reconstruction methods beyond those used for identify¬ 
ing the J/ijj or T. In particular, the decay H -7 J/ip'f 
may represent a viable probe of the Hcc coupling 0], 
which is sensitive to physics beyond the Standard Model 
(SM) BE], at the Large Hadron Collider (LHC). The 
expected SM branching fractions for these decays have 
been calculated to be B (H —> J/V’q) = (2.8 ± 0.2) x 
HT 6 , B[H^T{nS)i\ = (6.1±& 4 ,2.0±i;»,2.4±i;|) x 
10~ 10 [5]. No experimental information on these branch¬ 
ing fractions exists. These decays are a source of back¬ 
ground and potential control sample for the nonresonant 
decays H —» p + 7 . These nonresonant decays are sen¬ 
sitive to new physics [10] . 

Rare decay modes of the Z boson have attracted at¬ 
tention focused on establishing their sensitivity to new 
physics m- Several estimates of the SM branching frac¬ 
tion for the decay Z —» J/ijj 7 are available 021 - rH] with 
the most recent being (9.96 ±1.86) x 10 -8 [14]. Mea¬ 
suring these Z — y Q 7 branching fractions, benefiting 
from the larger production cross section relative to the 
Higgs case, would provide an important benchmark for 
the search and eventual observation of H —>• <2 7 decays. 
Additionally, experimental access to resonant Q 7 decay 
modes would also provide an invaluable tool for the more 
challenging measurement of inclusive associated Q 7 pro¬ 
duction, which has been suggested as a promising probe 
of the nature of quarkonium production in hadronic col¬ 
lisions 03 ns]. 


The decays Z —> Q 7 have not yet been observed, with 
the only experimental information arising from inclusive 
measurements, such as B (Z — > J/ipX) = ( 3 . 51 I 0 J 5 ) x 
10 ~ 3 and the 95% confidence level (CL) upper limits 
B[Z-+ T(nS)Y] < (4.4,13.9,9.4) x HT 5 , from LEP ex¬ 
periments 07H23- 

This Letter presents a search for decays of the recently 
observed Higgs boson and the Z boson to and 

T(nS)j final states. The decays J/-0 —► p + and 
T (nS) —)• p + are used to reconstruct the quarko¬ 
nium states. The search is performed with a sample 
of pp collision data corresponding to an integrated lu¬ 
minosity of 19.2 fb _1 (20.3 fb^ 1 ) for the J/V’T [T(nS) 7 ] 
analysis respectively, recorded at a center-of-mass energy 
y/s = 8 TeV with the ATLAS detector [22], described in 
detail in Ref. [25] . 

Higgs boson production is modeled using the 
POWHEG-BOX Monte Carlo (MC) event genera¬ 
tor [2IH25] . separately for the gluon fusion (ggF) and 
vector-boson fusion (VBF) processes calculated in Quan¬ 
tum Chromodynamics (QCD) up to next-to-leading or¬ 
der in a§- The Higgs boson transverse momentum (pt) 
distribution predicted for the ggF process is reweighted 
to match the calculations of Refs [22] [50], which in¬ 
clude QCD corrections up to next-to-next-to-leading or¬ 
der and QCD soft-gluon resummations up to next-to- 
next-to-leading logarithms. Quark mass effects in ggF 
production m are also accounted for. 

Physics beyond the SM that modifies the charm cou¬ 
pling can also change production dynamics and branch¬ 
ing fractions. In this analysis we assume the produc¬ 
tion rates and dynamics for a SM Higgs boson with 
mu = 125 GeV, obtained from Ref. [52], with an un¬ 
certainty on the dominant ggF production mode of 12%. 
The VBF signal model is appropriately scaled to account 
for the production of a Higgs boson in association with a 
W or Z boson or in association with a tt pair, correcting 
for the relative production rates and experimental accep¬ 
tances for these channels. Contributions from nonreso¬ 
nant H —> (Z*/ 7 *)q —> p + 7 decays are expected to be 
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negligible with respect to the present sensitivity [33M35] . 

The POWHEG-BOX MC event generator is also used 
to model Z boson production. The total cross section is 
estimated from Ref. [38], with an uncertainty of 4%. 

The Higgs and Z boson decays are simulated as a 
cascade of two-body decays. Effects of the helicity of 
the quarkonium states on the dimuon kinematics are ac¬ 
counted for in both cases. For Higgs and Z boson events 
generated using POWHEG-BOX, Pythia8.1 [37] [31] 
is used to simulate showering and hadronization while 
Photos [551 HD] is used to provide QED radiative correc¬ 
tions to the final state. The simulated events are passed 
through the full GEANT4 simulation of the ATLAS de¬ 
tector mm and processed with the same software used 
to reconstruct data events. 

The data used to perform the search in the J/rj) 7 chan¬ 
nel were collected using a trigger that required at least 
one muon with px > 18 GeV. The events used in the 
T ( nS ) 7 channel were collected with a trigger requiring 
an isolation muon with px > 24 GeV and a dimuon trig¬ 
ger with px thresholds of 18 and 8 GeV for each of the 
muons, respectively. Events are retained for analysis if 
they were collected under stable LHC beam conditions 
and the detector components were operating normally. 

Muons are reconstructed from inner-detector tracks 
combined with independent muon spectrometer tracks or 
track segments [55] and are required to have p^l > 3 GeV 
and pseudorapidity |p M | < 2.5. Candidate Q —► p + /i _ 
decays are reconstructed from pairs of oppositely charged 
muons consistent with originating from a common vertex. 
The liighest-px muon in a pair, called the leading muon in 
the following, is required to have pJ£ > 20 GeV. Dimuons 
with a mass, m^, within 0.2 GeV of the J/ip mass [T7] 
are identified as J/ijj —> /i + p“ candidates. In case both 
muons in the pair are within |p^| < 1.05, the said re¬ 
quirement is tightened to 0.15 GeV. Dimuons with 8.0 < 
< 12.0 GeV are considered as T (nS) —> p + p~ can¬ 
didates. The transverse momentum of each Q —> /i + p _ 
candidate, p 7 ,J , is required to exceed 36 GeV. 

Selected Q —> p + p~ candidates are subjected to isola¬ 
tion and vertex quality requirements. The sum of the px 
of the reconstructed inner-detector tracks and calorime¬ 
ter energy deposits within A R = \J (A</>) 2 + (A??) 2 = 0.2 
of the leading muon is required to be less than 10% of 
the muon’s px- The transverse momentum of the inner- 
detector track associated with the leading muon is sub¬ 
tracted from the sum and the subleading muon is also 
subtracted if it falls within the isolation cone. To reject 
backgrounds from 5-hadron decays, the measured trans¬ 
verse decay length L xy between the dimuon vertex and 
the primary pp vertex is required to be less than three 
times its uncertainty <jl x ■ In this case, the primary pp 
vertex is defined as the reconstructed vertex with the 
highest Y^iPri °f a h associated tracks used to form the 
vertex. 

Photon reconstruction is seeded by clusters of energy 


in the electromagnetic calorimeter. Clusters without 
matching tracks are classified as unconverted photon can¬ 
didates. Clusters matched to tracks consistent with the 
hypothesis of a photon conversion into an e + e _ pair are 
classified as converted photon candidates [Hj. Recon¬ 
structed photon candidates are required to have trans¬ 
verse momentum p^. > 36 GeV, pseudorapidity \rf\ < 
2.37, excluding the barrel/endcap calorimeter transition 
region 1.37 < |p 7 | < 1.52, and to satisfy the “tight” 
photon identification criteria |45| . To further suppress 
the contamination from jets, an isolation requirement is 
imposed. The sum of the transverse momentum of all 
tracks and calorimeter energy deposits within A R = 0.2 
of the photon direction, excluding those associated with 
the reconstructed photon, is required to be less than 8% 
of the photon’s transverse momentum. 

Combinations of a Q —> candidate and a photon, 

satisfying A^(p + /i _ ,7) > 0.5, are retained for further 
analysis. To improve the sensitivity of the search, the 
events are classified into four exclusive categories, based 
upon the pseudorapidity of the muons and the photon 
reconstruction classification. Events where both muons 
are within the region |ry M | < 1.05 and the photon is (is 
not) classified as a conversion constitute the “barrel con¬ 
verted” BC (“barrel unconverted” BU) category. Events 
where at least one of the muons is outside the region 
|p M | < 1.05 and the photon is (is not) classified as a con¬ 
version constitute the “endcap converted” EC (“endcap 
unconverted” EU) category. The number of candidates 
observed in each category following the complete event 
selection is shown in Table |T] 

The total signal efficiency (kinematic acceptance, trig¬ 
ger and reconstruction efficiencies) in the J/tjj 7 final 
state is 22% and 12% for the Higgs and Z boson de¬ 
cays, respectively. The corresponding efficiencies for the 
T(nS)7 final state are 28% and 15%. The m W7 reso¬ 
lution is similar for both the Higgs and Z boson decays 
and varies between 1.2% and 1.8%. The resolution 
is 1.4% and 2.4% for the barrel and endcap categories, 
respectively. 

The main source of background, referred to as inclusive 
QCD background, is dominated by inclusive quarkonium 
production where a jet in the event is reconstructed as 
a photon. For the T(n5)q final state, events contain¬ 
ing Z —> p + decays with final-state photon radiation 
(FSR) constitute a second source of background, a con¬ 
tribution which is found to be negligible in the J/ip 7 fi¬ 
nal state. The normalization of both of these background 
sources is extracted directly from a fit to data. The mod¬ 
eling of the inclusive QCD background shape, obtained 
with a data-driven approach, and of the Z —► pT back¬ 
ground shape, obtained from simulation, is described in 
the following two paragraphs. 

The background from inclusive QCD processes is mod¬ 
eled with a nonparametric data-driven approach using 
templates to describe the kinematic distributions. The 
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approach exploits a sample of loosely selected p, + p,~"/ 
events, around 2400 in the J/ij) 7 channel and around 
3200 in the Y(nS )7 channel. These control samples 
are formed from events satisfying the nominal Q 7 se¬ 
lection, but with relaxed dimuon and photon transverse 
momenta (p^ > 25 GeV and > 25 GeV) and isolation 
requirements (separate fractional calorimeter energy and 
track momentum isolation for the photon and dimuon 
system of less than 60%). Contamination of this sample 
from signal events is expected to be negligible. Prob¬ 
ability density functions (pdfs) used to model the pift*, 
Pt, Ap(p, + p~, 7 ) and Acf>(p + p~, 7 ) distributions of this 
control sample, independently for each category, are con¬ 
structed using Gaussian kernel density estimation [iHj . 
To account for kinematic correlations, the distributions 
of Ar](p + , 7 ) and A<j>{p + , 7 ) are estimated in 
eight exclusive regions of pjf 1 . In the case of the dimuon 
and photon isolation variables, correlations are accounted 
for by using two-dimensional histograms derived in five 
exclusive regions of p!£ M . The m M/1 distributions are mod¬ 
eled using Gaussian pdfs, with parameters derived from 
a fit to the control sample. In the T ( nS ) 7 channel, 
the data control sample is corrected for contamination 
from Z —>• p + p~ , y decays. The pdfs of these kinematic 
and isolation variables are sampled to generate an en¬ 
semble of pseudocandidates, each with a complete Q 7 
four-vector and an associated pair of correlated dimuon 
and photon isolation values. The nominal selection re¬ 
quirements are imposed on the ensemble and the surviv¬ 
ing pseudocandidates are used to construct templates for 
the kinematic distributions, notably the inclusive QCD 
background m A1M7 and piIjf distributions. 

The background from Z —> p + p .~7 decays is modeled 
with templates derived from a sample of simulated Z bo¬ 
son events with m in the T (nS) mass region. To vali¬ 
date this background model with data, the sidebands of 
the 771 ^^ 7 distribution in several validation regions, de¬ 
fined by relaxed kinematic or isolation requirements, are 
used to compare the prediction of the background model 
with the data. Good agreement within the statistical 
uncertainties is observed. 

The composition of the inclusive QCD background and 
the Z —> p + p “7 decay contribution is investigated with 
data. The details of the composition do not enter di¬ 
rectly the background estimation for this search, but the 
composition itself is a crucial input in feasibility studies 
for future searches or measurements, where projections 
of these backgrounds to different center-of-mass energies 
or luminosity conditions are needed. To facilitate this 
study, the selection requirements on m ^ and \L xy /<jL xy | 
are relaxed to include the sideband regions. In the J/ip') 
final state, a simultaneous unbinned maximum likelihood 
fit to the and | L xy /aL xy | distributions is performed. 
Once the simultaneous fit is performed, the composition 
of the subset of events satisfying the nominal and 
\L X y/&L xy | requirements is estimated. After the com¬ 


plete event selection, around 56% of the events originate 
from prompt J/%p production, 3% from non-prompt J/ip 
production (from 5-hadron decays) and 41% are combi¬ 
natoric backgrounds from nonresonant dimuon events. 

A separate simultaneous fit to the to /JM7 and dis¬ 
tributions of the same sample of candidate J/ip events 
finds no significant contribution from Z —> p + p ~7 de¬ 
cays, a conclusion that is also supported by a study based 
on simulated Z —> p + p.~ events. 

For the T(nS ')7 final state a simultaneous fit is per¬ 
formed to the m At/i7 and distributions. After the 
full event selection, inclusive T (nS) production accounts 
for 7% of events, 27% of the events are produced in 
Z —> / Lt + /r “7 decays, and 66 % of the events are asso¬ 
ciated with combinatoric backgrounds from nonresonant 
dimuon events. The contribution from Z —> p ~7 de¬ 
cays is in agreement with the MC expectation. 

TABLE I. The number of observed events in each analysis 
category. For comparison, the expected background yield is 
given in parentheses for the two m MM7 ranges of interest. The 
Higgs and Z boson contributions expected for branching frac¬ 
tion values of 10 -3 and 10 -6 , respectively, are also shown. For 
T(n 5 ) 7 , the 1S,'2S and 3S contributions are summed. 


S—1 

0 . 

Observed (Expected Background) 

Signal 

bO 

O 


Mass Range [GeV] 

Z 

H 

O 

All 

80-100 

115-135 

B [ 10 ~ 6 ] 

B [10~ 3 ] 

J/i> 7 

BU 

30 

9 (8.9A1.3) 

5 (5.0T0.9) 

1.29T0.07 

1.96±0.24 

BC 

29 

8 (6.0±0.7) 

3 (5.5±0.6) 

0.63±0.03 

1.06±0.13 

EU 

35 

8 (8.7A1.0) 

10 (5.8±0.8) 

1.37±0.07 

1.47±0.18 

EC 

23 

6 (5.6±0.7) 

2 (3.0T0.4) 

0.99A0.05 

0.93A0.12 

T(nS) 7 

BU 

93 

42 (39±6) 

16 (12.9±2.0) 

1.67T0.09 

2.6T0.3 

BC 

71 

32 (27.7±2.4) 

5 (9.7±1.2) 

0.79T0.04 

1.45±0.18 

EU 

125 49 (47±6) 

16 (17.8±2.4) 

2.24T0.12 

2.5T0.3 

EC 

85 

31 (31±5) 

18 (12.3±1.9) 

1.55T0.08 

1.60T0.20 


Trigger efficiencies and efficiencies for muon and pho¬ 
ton identification are determined from samples of Z —> ££, 
Z ££7 (£ = e,p) and J/ip —» p + p~ decays in 
data 03 03 . The systematic uncertainty on the ex¬ 
pected signal yield associated with the trigger efficiency 
is estimated to be 1.7%. The photon (both converted and 
unconverted) and muon reconstruction and identification 
efficiency uncertainties are estimated to be 0.5% (0.7%) 
and 0.4% (0.4%) for the Higgs boson {Z boson) signal, 
respectively. An uncertainty on the integrated luminosity 
of 2.8% is derived using the method described in Ref. j48l - 
The photon energy scale uncertainty, determined from 
Z —> e + e - and validated using Z —► ££7 decays 03 - 
is propagated through the simulated signal samples as a 
function of rf 1 and p'/. The uncertainty associated with 
the description of the photon energy scale in the simu¬ 
lation is found to be less than 0 . 2 % of the three-body 
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invariant mass while the uncertainty associated with the 
photon energy resolution is found to be negligible rela¬ 
tive to the overall three-body invariant mass resolution. 
Similarly, the systematic uncertainty associated with the 
muon momentum measurement is determined using data 
samples of J/ip — > p + p~ and Z -+ p + p~ decays and 
validated using T (nS) -+ p + p~ decays pf5] . For the 
Pt range relevant to this analysis, the systematic uncer¬ 
tainties associated with the muon momentum scale are 
negligible. 

The uncertainty in the shape of the inclusive QCD 
background is estimated through the study of variations 
in the background modeling procedure. The shape of 
the pdf is allowed to vary around the nominal shape 
within an envelope associated with shifts in the pif/ 1 and 
p'ry distributions. Furthermore, a separate background 


model, generated without removing the contamination 
from Z -+ p + decays, provides an upper bound on 
potential mismodeling associated with this process. 

Results are extracted by means of a simultaneous 
unbinned maximum likelihood fit, performed to the 
selected events with 30 GeV < m M/i7 < 230 GeV 
separately in each of the analysis categories. In the 
final state, the fit is performed on the m W7 and 
p distributions, while for the T(nS)"f candidates 
a similar fit is performed using the m^-y, p ^ 1 , and 
distributions. The latter distribution provides 
discrimination between the three T(nS) states and 
constrains the Z —> p + p~ 7 background normalization. 
No significant Z —> Qj or H —> Q 7 signals are observed, 
as shown in Figs, [l] and [ 2 ] 




FIG. 1. (color online) The and p 7 M7 distributions of the selected J/ipl candidates, along with the results of the unbinned 

maximum likelihood fit to the signal and background model (S+B fit). The error bars on the data points correspond to the 
statistical uncertainties. The Higgs and Z boson contributions as expected for branching fraction values of 10~ li and 1CF 6 , 
respectively, are also shown. 
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FIG. 2. (color online) The m^-y, p 7 M7 and m distributions of the selected T (nS) 7 candidates, along with the results of the 
unbinned maximum likelihood fit to the signal and background model (S+B fit). The error bars on the data points correspond 
to the statistical uncertainties. The Higgs and Z boson contributions as expected for branching fraction values of 10~ J and 
10~ 6 , respectively, for each of the T (nS) are also shown. 
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Upper limits on the branching fractions for the Higgs 
and Z boson decays to J/-i />7 and T(nS)^ are set us¬ 
ing the CL S modified frequentist formalism [50] with the 
profile likelihood ratio test statistic [51] . The expected 
SM production cross sections are assumed for the Higgs 
and Z bosons. The results are summarized in Table im 
The 95% CL upper limit on the branching fraction 
for H —> J/ip'l decays corresponds to about 540 times 
the expected SM branching fraction. The upper limits 
on the Z —> J/ip') and Z —> Y(nS)~f branching frac¬ 
tions significantly constrain the allowed range of values 
obtained from theoretical calculations Upper 

limits are also set on the combined branching fractions 
B[H-> T(nS) 7 ] < 2.0 x 10 " 3 and B [Z -> T(nSb] < 
7.9 x 1CU 6 , where the relative contribution of each final 
state to the potential signal is profiled (allowed to float 
to the values that maximize the likelihood) during the 
fit. 

TABLE II. Expected and observed branching fraction limits 
at 95% CL for yfs = 8 TeV. The ±lcr fluctuations of the 
expected limits are also given. For the Higgs decay search, 
limits are also set on the cross section times branching fraction 
a (pp -4 H) x B (H -4 Qy). 



95% CL Upper Limits 



T(1S) 

T(2 S) 

Y(3S) 

E n ?(nS) 

B(Z^ Q 7 ) [ 10- 6 ] 

Expected 

2 0 +1 ° 
Z>U -0.6 

4.9 

6-2±f;l 

5-4 t 2 d 

8-8 

Observed 

2.6 

3.4 

6.5 

5.4 

7.9 

B(H -4 Q 7 ) [ 10 -3 ] 

Expected 

1 9 +O .6 
—0.3 

1 8 +0,9 
1,O -0.5 

2 1 

Z *' L — 0.6 

i-s±S:S 

O c + 1-3 
z -°-0.7 

Observed 

1.5 

1.3 

1.9 

1.3 

2.0 


° (pp 

-4 H) x B (H -4 Q 7 ) [fb] 


Expected 

26±]. 2 

38±J? 

45i?^ 

38±« 

54 til 

Observed 

33 

29 

41 

28 

44 


In conclusion, the first search for the decays of the SM 
Higgs and Z bosons to J/xj) 7 and Y(nS) 7 (n = 1,2,3) 
has been performed with y/s = 8 TeV pp collision data 
samples corresponding to integrated luminosities of up to 
20.3 fb -1 collected with the ATLAS detector at the LHC. 
No significant excess of events is observed above the back¬ 
ground. In the J/'0 7 final state the 95% CL upper limits 
on the relevant branching fractions for the SM Higgs and 
Z bosons are 1.5 x 1CP 3 and 2.6 x 10 —6 , respectively. The 
corresponding upper limits in the T(1S, 2S, 3S) 7 chan¬ 
nels are (1.3,1.9,1.3) x 10 ~ 3 and (3.4, 6.5, 5.4) x 10~ 6 , for 
the SM Higgs and Z bosons respectively. These are the 
first experimental bounds on exclusive Higgs and Z bo¬ 
son decays to final states involving quarkonia. 
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